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Abstract. We demonstrate that in a certain gauge the Lax matrices of the 
rational and hyperbohc Ruijsenaars-Schneider models have a quadratic r-matrix 
Poisson bracket which is an exact quadratization of the linear r-matrix Poisson 
bracket of the Calogero-Moser models. This phenomenon is explained by a geo- 
metric derivation of Lax equations for arbitrary flows of both hierarchies, which 
turn out to be governed by the same dynamical i?-operator. 



1 Introduction 



In the recent years the interest in the Calogero-Moser type of models 
is considerably revitalized. One of the directions of this recent development 
was connected with the notion of the dynamical r-matrices and their inter- 
pretation in terms of Hamiltonian reduction 0-0]. Very recently, this line 
of research touched also the so-called Ruijsenaars-Schneider models |T3|,[p!4 



which may be seen as relativistic generalizations of the Calogero-Moser ones 



In the present paper the same subject as in [IJ] is handled, namely, 
the rational and hyperbolic Ruijsenaars-Schneider models. However, the 
results are somewhat different, and, as we hope, somewhat more beautiful. 
The difference is due to another gauge we choose for the Lax matrix of the 
Ruijsenaars-Schneider models. A striking circumstance comes out when us- 
ing our gauge, namely that the both classes of models are governed by one 
and the same dynamical i?-operator. This seems to pass unnoticed in the 
existing literature and is explained in Sect. 4. The computations presented 
there are hardly new, at least for the non-relativistic case, but we could 
not find in the literature the main message following from these computa- 
tions, namely that they give Lax reprentations for an arbitrary flow of the 
corresponding hierarchies, and hence are perfectly suited for guessing (not 
proving!) the correct r-matrix ansatz. Sect. 3 contains the main result of the 
paper, namely a quadratic r-matrix Poisson structure for the rational and 
hyperbolic Ruijsenaars — Schneider models. We compare our results with the 
previously known ones in the Sect. 5. Besides, for the convenience of a general 
reader we give a short review of relevant notions from the r-matrix theory 
in the Sect. 2. Sect. 6 is devoted to some problems arising from our results. 



2 General framework 

In this section we recall some basic notions about integrable hierarchies and 
their r-matrix theory. Our formulations result from the observations on the 
large "experimental material" collected in the last decades of research in this 
area. 

Let V he a. Poisson manifold; in fact we consider here only the simplest 
case of a symplectic space R^^{x,p} with canonically conjugated coordinates 
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y and p = (pi, . . . ,Pn)'^, so that the Poisson brackets of the 
coordinate functions read: 



{xk,Xj} = {pk,Pj} = 0, {xk,Pj} = 5kj. (2.1) 

Let H{x,p) be a completely integrable Hamiltonian, i.e. suppose that the 
hamiltonian system 

x = {x,H}, p={p,H} (2.2) 

posesses functionally independent integrals in involution. Then (|2^ ) usu- 
ally (probably always) admits a Lax representation, i.e. there exist two maps 
L : V ^ g and M : V ^ g into some Lie algebra g such that (|2.2|) is 
equivalent to 

L=[M,L]. (2.3) 

In the cases we are dealing with in this paper, g = gl{N), and it carries some 
additional structures. In particular, it is an assosiative algebra with respect to 
the usual matrix multiplication, and admits a non-degenerate scalar product 
{U, V) = tY{UV). In other cases g could be a more complicated algebra, for 
example a loop algebra; then one would speak about Lax representation with 
a spectral parameter. 

An important observation is that usually 

H{x,p) = ip{L), 

where (p{L) is an Ad-invariant function on g. This is related to the fact that 
integarble systems appear not separately, but are organized in hierarchies. 
Namely, to every invariant function <f{L) there corresponds a Lax equation 
of the form ( p.3| ). Moreover, there often holds the following relation: 

M = R{V^{L)), (2.4) 

where R : g g is a linear operator, depending, may be, on some of the 
coordinates on V (it is called then dynamical) . We shall call R an R-operator 
governing the corresponding hierarchy. (Recall that the gradient Vip{L) of a 
smooth function ip on an algebra g with a scalar product {U, V) is defined 
by the relation 



(Vv.(L),f/)= ^^(L + ef/) 



e=0 
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and for an Ad-invahant (p one has [Vv9(L), L] =0). 

Sometimes the relation (2^) has actually another form: 



M = R{LVip{L)). (2.5) 



Of course, in the general setting, where R in ( p.4|) can be dynamic, the 
equation (|2.5| ) can be seen as M = lZ{Vip{L)) with an operator Ti{-) = 
R{L-). However, if, as it is often the case in applications, an operator R 
in (|2.5| ) has some special properties (e.g. is independent on some or all 
dynamical variables), then it is advantageous to consider this particular case 
on its own rights. We shall call R in ( |2.5|) also an R-operator governing the 
corresponding hierarchy, keeping in mind differnce between ( |2.4| ) and ( |2.5| ). 



An r-matrix theory |T^, provides a sort of explanation of relations 
(|2.5| ). Namely, the formula (|2.4| ) is usually a consequence of a more 
fundamental fact, namely that L = L{x,p) form a Poisson submanifold in g 
equipped with a so called linear r -matrix bracket. This is expressed in the 
formula 

{Lf L} = [J®L,r] - [L® J,r*]. (2.6) 
The A^^ X A^^ matrix r corresponding to the operator R is defined as 

N 

i,j,k,m=l 

where 

rij^km = {R{Eji),Emk) = coeS. by E^m in R{Eji). (2.7) 
and r* corresponds in the same way to the operator R*, so that 

AT 

i,j,k,m=l 

In case of constant (non-dynamical) -R-operators a sufficient condition for 
( |2.6| ) to define a Poisson bracket is given by the so called modified Yang- 
Baxter equation In case of dynamical /^-operators the corresponding 
theory is less developed; the most general known sufficient conditions are 
given in [||. The lax representation ( p.3|) with M given in (|2.4| ) is a conse- 
quence of of (|2.6| ) for an arbitrary Arf-invariant function <f{L), 
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Analogously, the formula ( p.5| ) is usually a consequence of a more funda- 
mental fact that L = L{x,p) form a Poisson submanifold in g equipped with 
a so called quadratic r -matrix bracket. The most general quadratic bracket 
is given by 

{LfL} = {L®L)ai-a2{L®L) + {I®L)si{L®I)-{L®I)s2{I®L). (2.8) 
Such general quadratic r-matrix structures were discovered several times in- 



dependently [IT3], |1T8[, [0; see ||T9[ for an application to the Toda and rela- 
tivistic Toda hierarchy. 

In (|2.8| ) the matrices ai, 02, Si, S2 correspond in a canonical way to some 
linear (in principle, dynamical) operators Ai, A2,Si, S2 and satisfy conditions 

a* = — fli, 02 = — a2, 3*2 = Si, (2-9) 

and 

ai + si = 02 + S2 = r. (2.10) 

The first of these conditions assures the skew-symmetry of the Poisson 
bracket (|2.8|) , and the second one garantees that the Hamiltonian flows with 
invariant Hamiltonian functions ^{L) have the Lax form ( p.3| ) with the M- 
matrix (|2.5|) . If (|2.10|) is satisfied, we call the bracket (|2.8| ) a quadratization 
of the bracket ( |2.6| ). In the case of constant operators a sufficient condition 
for ( p.8| ) with ( p.9|) , ( p. 10 ) to be a Poisson bracket is validity of the modified 



Yang-Baxter equation for three operators i?, Ai, A2; nothing is known for 
dynamical case. 

It should be remarked that, while in the linear case the correspondence 
between operator R in ( |2.4| ) and matrix r in ( p.6| ) is rather unambigous, 
the situation is quite different in the quadratic case. There exist in principle 
infinitely many possibilities of choice of ai, 02, si, $2 in (p78|), satisfying ( p.9| ), 
( |2.10|) and resulting in the same operator R in (|2.5| ). All such quadratizations 



are parametrized by one skew-symmetric matrix, for example, by ai, because 
a2 = r — r* — ai, Si = r — ai, S2 = r* + ai. Hence finding a quadratic r-matrix 
Poisson structure for a given Lax matrix L is a non-trivial entertainment even 
if the i?-operator governing the whole hierarchy is known. Guessing a correct 
quadratization remains more or less a matter of art. The present paper 
is devoted to finding such quadratic structure for rational and hyperbolic 
Ruisenaars-Schneider models. 
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3 Ruijsenaars— Schneider models 
and their r— matrix formulation 

The hyperbolic relativistic Ruijsenaars-Schneider (RS) hierarchy is described 
in terms of the Lax matrix 

L = Lns{x,p) = E ■ .r^^^\ , hEkr (3.1) 

smh(xfc - Xj + 7) 

Here 7 is a parameter of the model, usually supposed to be pure imaginary. 
We use an abbreviation 

h = exp(p.) n f 1 - -#^^1 ' (3-2) 

j^k \ smh {xk-Xj)J 

so that in the variables (x, b) the canonical Poisson brackets ( p.l| ) take the 
form 

{xk, Xj} = 0, {xk, bj} = bkSkj, {bk, bj} = nkjbkbj (3.3) 

with 

TTkj = coth(a;j - + 7) - coth(a;fc - xj + 7) + 2(1 - 6kj)coth.{xk - xj). (3.4) 

The Hamiltonian of the RS model proper (i.e. of the simplest member of 
the hierarchy) is given by 

N 

H{x,p) = E 6fc = tr L{x,p). 

k=l 

This model admits a so called non-relativistic limit, achieved by rescaling 
p I— > Pp, 7 I— >• /57 and subsequent sending P 0. In this limit we have 
-^RS = I + PLcM + 0(/5^), where the Lax matrix of the rational Calogero- 
Moser (CM) hierarchy is introduced: 

L = Lcm{x,p) = 22PkEkk + y2^-rf \Ekj- (3.5) 

^. smh(xfc-Xj) 

The Hamiltonian of the CM model proper is given by 

1^1 7^ 1 
H(x, p) = - H — n- = -tr L^(x, p). 

^ k=l ^ k+3 ^1^"^ -^3 1 ^ 
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A simple, but important particular case of these models constitute the 
rational ones, which are obtained by rescaling x i— >■ fix, 7 //7 and sending 
/X — > 0. In this limit one gets the Lax matrix of the RS hierarchy: 

^ 7 

L = Lns{x,p) = V ■ — bjEkj, (3.6) 

where 

h = expipk) n f 1 - r-^] ■ (3.7) 
The canonical Poisson brackets in terms of (x, b) are still given by (|3.3|) with 

L^ + !(i^, (3.8) 

Xj - Xk + J Xk- Xj +-f Xk - Xj 

As a Lax matrix of the CM hierarchy one gets 

L = Lcm{x,p) = Y,PkEkk + Yl ^kj- (3.9) 

Now we are in a position to formulate the main result of this paper. 
Theorem. For the Lax matrices of the relativistic models ( |3.6| ), (|3.1| ) 
there holds a quadratic r -matrix ansatz (|2.8| ) wi/i t/ie matrices 

ai = a + w, Si = s — w, 

a2 = a + s — s*— w, S2 = s* + w , 
where in the rational case 

a = Yl Ejk^Ekj, (3.10) 

k-tj Xk Xj 

s = -J2 — ^Ejk^Ekk, (3.11) 

k^j Xk Xj 

w = J2 -^—Ekk ® Ej„ (3.12) 

k^j ■'^k Xj 



6 



and in the hyperbolic case 

a = ^ coth(xfc - Xj)Ejk (S) Ekj, (3.13) 

^ = -E-i:7-^ ^E,k®Ekk, (3.14) 

smh(xfc - Xj) 

W7 = ^ COth(xA: — Xj)Ekk <S> Ejj. (3.15) 

Let us stress three remarkable properties of the found r-matrix structure. 

• All matrices a, s, w are dynamical, but depend only on the coordinates 
Xk, not on the momenta pk- 

• All matrices a, s, w do not depend on the parameter 7 of the model. 

• The most important and striking fact is that the structure found is an 
exact quadrat iz at ion of a linear r-matrix bracket with 

r = a + s, 

which is just the r-matrix of the non-relativistic CM model found 
by Avan-Talon in 0. So, the rational and hyperbolic RS hierarchies 
turn out to be governed by the same i?-operators as the rational and 
hyperbolic CM hierarchies, respectively. 

The proof of the Theorem above is a matter of rather direct computa- 
tions, and is therefore omitted. The next section is devoted to a geometric 
derivation of the /^-operators for the RS and CM hierarchies. This deriva- 
tion is independent on the our Theorem, and provides some explanation of 
the third property mentioned above. 



4 Derivation of i?— operator for CM and RS 
hierarchies 

An important tool to investigate the hyperbolic CM and RS models is, along 
with the Lax matrices, an auxiliary diagonal matrix 

N 

X = diag(exp(a;i), . . . , exp(xAr)) = ^ exp{xk)Ekk- (4.1) 

k=l 
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The fundamental commutation relation, describing the structure of the Lax 
matrix L in terms of a given diagonal matrix X, reads for the CM model: 

XLX-^ - X-^LX = 2-fJ2 ^kj = 27(ee'^ - /), (4.2) 

and for the RS model: 

e''XLX-^ - e-^'X-^LX = 2sinh(7)eb'^. (4.3) 
Here / stands for the N x N unity matrix, e = (1, . . . , 1)"^, and b = 

(6i,...,M^- 

We shall derive the i?-operators of these models from the results about 
their explicit solution obtainable from |1|], 0. These results may be formu- 
lated as follows. Let y:>{L) be an Ad-invariant function on gl{N), and take 
its value on one of the Lax matrices Lqm, -Z^rs as a Hamiltonian function 
for the corresponding model. Then the quantities exp(2a;fc(t)) are just the 
eigenvalues of the matrix 

X2exp(2t/(Lo)), 

where 

f{L) = Vip{L) for the non — relativistic CM case, (4-4) 

f{L) = LS/ip{L) = Vip{L)L for the relativistic RS case. (4.5) 

/^From this statement the Lax equations for the corresponding flows can 
be derived. Indeed, define the evolution of the matrices X, L by the equations 

= X\t) = VXl exp{2tf{Lo))V-\ (4.6) 

L = L{t) = VLoV-\ (4.7) 

Let us explain, in which sence these equations define an evolution. The 
matrix X"^ = X'^{t) consists of eigenvalues of the matrix Xq exp(2t/(Lo)), and 
the matrix V = V{t) is a diagonalizing one. (It is easy to see that the matrix 
Xq exp(2t/(Lo)) is similar to a self-adjoint one, if 7 is pure imaginary). If we 
fix the ordering of Xk (for example, Xi < . . . < x^), then the only freedom in 
the definition of is a left multiplication by a diagonal matrix. We fix now 
V by the condition 

VXoe = Xe, (4.8) 
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and show that this assures that the corresponding requirements ( |4.2| ) and 
(^4.3|) are satisfied for all t, provided they were satisfied for t = 0. 
Indeed, we have for the CM case: 

XLX-^ - X-^LX = X-^VXo [XqLoX^^ - X^^LoXq) X^^V'^X 

= 27X-VXo(ee^ - I)Xo^V-^X = 27 (X" VXcee^Xg^V^^X - /) . 

Since the diagonal entries of the matrix on the left-hand side vanish, we see 
that (J^ implies e^X^^V-^X = e^, which proves the commutation relation 
(PD for all t. 

Analogously, for the RS case we have: 

e^XLX-^ - e~^X-^LX 

= X-'VXo {e^XoLoX^' - e~^Xo-^LoXo) X,'V-'X 
= 2sinh(7)X~VXoeb^Xo-V-^X. 

so that denoting bQXQ"^V"~^X = b"^, we see that ( [4.8| ) implies the validity of 
the commutation relation (O) for all t. 

/^From this point the derivation of the evolution equations for L, X is 
identical for both the non-relativistic and the relativistic cases. Differentiat- 
ing Q, (|;§, we get: 

L=[M,L], (4.9) 
(so that the evolution of L is governed by a Lax equation), 

2XX = [M, X^] + 2X^f{L), (4.10) 

where 

M = VV-\ (4.11) 

In order to find the matrix M explicitly, consider first the off-diagonal part 
of (^4.101) , which implies: 

= exp(3;fc-x,) ^ _ ))^(^) ^ ^ ^. (412) 

smh(Xfc — Xj) 

The normalizing condition following from ( |4.8| ), ( |4.11| ) reads: 

MXe = Xe, (4.13) 
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hence 

Mkk = Xk-Yl Mkj exp(xj - Xfc). 

An expression for Xk can be read off the diagonal part of ( [4.10 ): 

Xk = f{L)kk- 

Substituting this in the previous formula and using (|4.12| ), we get: 

Mkk = f{L)kk - E ■ 3^^"' y (4.14) 
^^smh[xk-Xj) 

Finally, notice that we can redefine M — f{L) as a new M (this does not 
influence the equations of motion described by the Lax pair), which results 
in more convenient expressions 

Mkj = coth(xfc - Xj)f{L)kj, k 7^ j, (4.15) 

M.. = -E-#^- (4.16) 
smh(a;fc - xj) 

Analogous results for the rational case may be obtained by the limiting 
process, or derived in parallel. An auxiliary diagonal matrix is then given by 

N 

X = diag{xi,...,XN) = J2^kEkk- (4.17) 

k=l 

Given the diagonal matrix X, the structure of the Lax matrix L is completely 

described by the following fundamental commutation relations: for the CM 
case 

XL - LX = 7 E Ekj = 7(ee^ - /), (4.18) 

and for the RS case 

XL- LX + -fL = 7eb^. (4.19) 

The results about an exphcit solution of these models |]1|, [0] read: let 
(/9(L) be an Ad-invariant function on gl{N), and take its value on one of 
the Lax matrices ( |3.9| ), ( |3.(j| ) as a Hamiltonian function of the corresponding 
model. Then the quantities Xk{t) are just the eigenvalues of the matrix 

Xo + tf{Lo), 
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where ( [4.4| ), ( |4.5| ) still hold. Defining the evolution of the matrices X, L by 
the equations 

X = X(t) = V{Xo + tf{Lo))V~\ (4.20) 

L = L{t) = VLoV-\ (4.21) 

we see that in order to assure that the commutation relations (|4.18| ) and 
( ^.19| ) are preserved in the dynamics, one has to normalize the matrix V by 
the condition 

Ve = e, (4.22) 

Now the calculations are again identical for both the CM and the RS 
cases. Differentiating ( [4.21| ), ( [4.20| ), we get: 

L=[M,L], (4.23) 

(so that the evolution of L is governed by a Lax equation), 

X=[M,X] + f{L), (4.24) 

where 

M = VV-\ (4.25) 

Differentiating ([4.22|) gives: 



Me = 0, 



(4.26) 



Now we can find the matrix M explicitly. From ( |4.24| ) we immediately 
obtain the off-diagonal entries of the matrix M: 



f{L)kj 



The normalizing condition ( f4.26| ) implies: 



Mi 



kk 



' E Mkj 



f{L)k, 



Xk 



(4.27) 



(4.28) 



According to (|4.4| ), ( [4.5|) , we see that the following statement is proved. 

Proposition. The general flows of the (rational or hyperbolic) CM and 
RS hierarchies with Ad-invariant Hamiltonians tf{L) have Lax representa- 
tions i\2.Si ) with the M -matrices given by i\2.4 ), (^A ), respectively. The R 
operator is one and the same for the CM and RS cases and is given by: 



R = A + S, 



(4.29) 
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where A is a skew-symmetric operator on gl{N), and S is a non-skew- 
symmetric one, whose image consists of diagonal matrices. For the rational 
models 

A{Ekj) = ^—^Euj, (4.30) 

X]^ Xj 

S{Ekj) = -^—^Ekk, (4.31) 
and for the hyperbolic models 

A{Ekj) = (1 - 5kj)coth{xk - Xj)Ekj, (4.32) 

S{E,,) = - . ^ E,,. (4.33) 

smh[Xk — Xj] 



Obviously, operators A, S canonically correspond to the matrices a, s 
from (1^3), (|T|) and (|3T3|) , (CT ), respectively. 



5 Comparison with the previous results 

As pointed out in the Introduction, the r-matrices for the RS models were 



previuosly discussed by Babelon-Bernard in |13| and by Avan-RoUet in [|I4 
In both papers another gauge for the Lax matrix is chosen, namely a self- 
adjoint one: 

L=j: ■ A hh.Y'^E,, (5.1) 

smh(a;fc - Xj + 7) 

Let us reformulate our results for this gauge. Calculations analogous to those 
presented in Sect. 4 show that the M-matrix for the flow of RS hierarchy with 
a Hamiltonian function ^{L) is given in the gauge ( |5.1|) by 

Mkj = coth(a;fe - Xj)f{L)kj, j ^ k, 

M.. = -li:fM"'Z(%-f(%. (5.2) 

2^^\hJ smh(xfc-Xj) 

where f{L) = LVcp^L), as in ([4.5|). Accordingly, the matrices ai, 02, si, 
S2 in a quadratic Poisson bracket (|2.8| ) for the matrix ( |5.1|) depend with 
necessity not only on Xk but also on the momenta pk- After direct but 
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tedious calculations one can get a gauged version of the bracket from our 
Theorem in the form fE.SI) with 



ai = a + ^{ui-ul) +v, Si = ^{u2 + ul) -V, 

1 1 

a2 = a + -{u2 - U2) - V, S2 = -{ui + ul) + V, (5.3) 

where 

Ui=Y^ UkjEkj ® Ekk, U2 = -Y^ UkjEjk ® Ekk, V = Y1 '"kjEkk ® Ejj, 
the coefficients of these matrices being given by 

Vkj = ^(coth(a;fc — Xj + 7) — coth(xj — + 7) + 2(1 — 6jk)coth.{xk — Xj 



and 



Uki 



\^k/ smh.{xk — Xj) 

Let us stress that these matrices, first, depend on momenta pk, second, de- 
pend on the parameter 7, and third, do not reproduce the i?-operator govern- 
ing the CM hierarchy. This demonsrates how crucially depend the properties 
of r-matrices on the choice of gauge for Lax matrix. 

5.1 On the Avan— Rollet r— matrix 

Avan-Rollet found a linear r-matrix structure of the type ( p.6| ) for the Lax 



matrix (|5.1| ) . They looked for an r-matrix with entries being linear combina- 
tions of elements Lkj of the Lax matrix ( |5.1| ) with coefficients depending only 
on coordinates Xk, but not on the momenta pk- According to the remarks 
above, such linear combinations can not be cast in the form 

r = (L(g)/)ri + r2(L(g)/), (5.4) 

necessary for putting linear r-matrix structure ( |2.6|) into quadratic form 



( |2.8|) . Indeed, in order to have ( [5.4| ) one is forced to admit coefficients de- 
pending on momenta. 
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Actually, after some manipulations based on the explicit expressions of 
Lkj, the Avan-RoUet r-matrix can be demonstrated to be equal to ( ^.41) with 

n = -{a + ui + v), r2 = -{a + U2-v), 

which gives back the expressions (|5.3| ). Hence our results and those by Avan- 
RoUet are in a sort of "hidden" equivalence. 

5.2 On the Babelon— Bernard r— matrix 

Babelon-Bernard consider the hyperbolic RS system with a specific value of 
parameter 7 = i7r/2 (which, by the way, precludes the possibility of both 
the non-relativistic and the rational reductions). They use the gauge (|5.1| ) 
and obtain a quadratic Poisson bracket with ai = 02 = |a, si = S2 = \a7 , 
where r denotes the transposition in the first factor of tensor products. In 
particular, all r-matrix objects depend only on x^'s. This seems to be in 
contradiction with the remark after ( |5.2|) . However, a closer look at (|5.2| ) 
solves this paradox. For 7 = m /2 the Lax matrix L ( |5.1| ) is symmetric, which 
forces diagonal entries Mkk of the matrix M to vanish. It means that for this 
particular case the i?-operator is equal symply to A, or, taking into account 
the symmetry of L, to the + AoT), where T stands for the transposition 
operator. The structure found by Babelon-Bernard is the simplest possible 
quadratization of a twisted linear r-matrix bracket with r = |(a + a^). So, 
it essentially takes care of specific properties of the model for 7 = and 
is not covered by our general construction. 



6 Conclusions 

The results of this paper surely constitute only one link in a long chain of (al- 
ready achieved and still hypothetical) results concerning the RS models. The 
whole work made for their non-relativistic counterparts should be repeated, 
or, better, generalized. 

• After this paper was submitted for publication, two apparently differ- 
ent (spectral parameter dependent) quadratic r-matrix structures were 



found for the elliptic RS hierarchy |2B|]. The structure found in the 
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second of these preprints serves as a direct generalization of the present 
results: it turns out that the elliptic RS hierarchy is also governed by 
the same /^-operator as the elliptic CM one. 



The spin Ruijsenaars-Schneider models, introduced recently in 



should be also put into the r-matrix framework, as it was done for the 
spin Calogero-Moser systems in |p. 

Another important point is a derivation of our Theorem in the frame- 



work of Hamiltonian reduction; to this end the results in ||TI[ should 
be used and further developed. 

• It would be rather important to investigate the dynamical analogs of the 
modified Yang-Baxter equation assuring the Poisson bracket properties 
of the ansatz (|2.8|) . The investigations of these objects were started in 
for the linear ansatz ( p.6| ) and are expected to unveil new interesting 
structures in the case of quadratic brackets. 

Further, our Theorem being established, situation with the Calogero- 
Moser type models becomes perfectly analogous to the situation with the 
Toda-like ones. Namely, the transition from the non-relativistic to the rela- 
tivistic systems corresponds to the transition from a linear r-matrix Poisson 
structure to its quadratisation. (See [^, [|1^] for the Toda case). A deeper 



understanding of this phenomenon is highly desirable. 
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